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INTRODUCTION 



In the chiral hmit the QCD Lagrangian has an important global chiral symme- 
try: Specifically, for QCD with three colors and Nf massless flavors, the invariance 
group acting on the quark flelds is SU{Nf) x SU{Nf) x [/^(l). In spite of the light- 
en I ness of the up and down quark, no signatures of this symmetry have been observed in 
nature. On the contrary, the fact that the pion mass is much lighter than the other 
hadron masses indicates that SU{Nf) x SU{Nf) is spontaneous broken with the pi- 
ous as Goldstone bosons. The breaking of Ua{^) follows from the absence of parity 



D . doublets. For example, the mass of the 6 is very different from the mass of the pion. 

^ ' The description of the broken phase requires the introduction of an order parameter 

which is zero in the restored phase. For the chiral phase transition the order parameter 
^ ■ is the chiral condensate (ipip)- Large scale lattice QCD simulations performed during 

■ the past decade show that for two light flavors the value of the chiral condensate is 

approximately —{220MeV)^ at zero temperature and vanishes above a critical tem- 
perature of about 140 MeV (see reviews by DeTar^, Ukawa^ and Smilga^ for recent 
results on this topic). As was in particular pointed out by Shuryak"^, not necessarily 
both symmetries are restored at the same temperature^. This may lead to interesting 
observable consequences. 

According the Banks-Casher formula^, the chiral condensate is directly related to 
the spectrum of the Euclidean Dirac operator near zero virtuality (here and below, we 
use virtuality for the value of the Euclidean Dirac eigenvalues). However, the eigenval- 
ues fluctuate about their average position as the gauge flelds vary over the ensemble 
of gauge fleld conflgurations. The main question that will be addressed in these two 
lectures is to what extent such fluctuations are universal. Inspired by the study of 
spectra of complex systems'', we will conjecture that spectral fluctuations on the scale 
of a typical eigenvalue spacing are universal^. If this conjecture is true, such spectral 
fluctuations can be obtained from a broad class of theories with global symmetries of 
QCD as common ingredient. We will derive them from the simplest models in this 
class: chiral Random Matrix Theory (chRMT). 

A flrst argument in favor of universality in Dirac spectra came from the analysis 
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of the finite volume QCD partition function^. As has been shown by Gasser and 
Leutwyler^°, for box size L in the range 

l/A<L<l/m^, (1) 

(A is a typical hadronic scale and m^r is the pion mass) the mass dependence of the 
QCD partition function is completely determined by its global symmetries. As a con- 
sequence, fluctuations of Dirac eigenvalues near zero virtuality are constrained by, but 
not determined by, an infinite family of sum rules^ (also called Leutwylcr-Smilga sum 
rules). For example, the simplest Lcutwylcr-Smilga sum rule can be obtained from the 
microscopic spectral density^ (the spectral density near zero virtuality on a scale of a 
typical eigenvalue spacing). On the other hand, the infinite family of Leutwyler-Smilga 
sum rules is not sufficient to determine the microscopic spectral density. The addi- 
tional ingredient is universality. A priori there is no reason that fiuctuations of Dirac 
eigenvalues are in the same universality class as chRMT. Whether or not QCD is inside 
this class is a dynamical question that can only be answered by full scale lattice QCD 
simulations. However, the confidence in an affirmative answer to this question has been 
greatly enhanced by universality studies within chiral Random Matrix Theory. The aim 
of such studies is to show that spectral fluctuations do not depend on the details of 
the probability distribution. Recently, it has been shown that the microscopic spectral 
density is universal for a wide class of probability distributions^^' We will 

give an extensive review of these important new results. 

Because of the C/a(1) symmetry of the Dirac operator two types of spectral fiuc- 
tuations can be distinguished. Spectral fluctuations near zero virtuality and spectral 
fluctuations in the bulk of the spectrum. The fluctuations of Dirac eigenvalues near 
zero virtuality are directly related to the approach to the thermodynamic limit of the 
chiral condensate. In particular, knowledge of the microscopic spectral density provides 
us with a quantitative explanation^^ of finite size corrections to the valence quark mass 
of dependence of the chiral condensate^^. 

Recently, it has become possible to obtain all eigenvalues of the lattice QCD 
Dirac operator on reasonably large lattices^^' making a direct verification of the 
above conjecture possible. This is one of the main objectives of these lectures. This is 
easiest for correlations in the bulk of the spectrum. Under the assumption of spectral 
ergodicity^^ , eigenvalue correlations can be studied by spectral averaging instead of 
ensemble averaging^^' On the other hand, in order to study the microscopic spectral 
density, a very large number of independent gauge fleld conflgurations is required. First 
lattice results conflrming the universality of the microscopic spectral density have been 
obtained recently^". 

At this point I wish to stress that there are two different types of applications of 
Random Matrix Theory. In the flrst type, fluctuations of an observable are related to 
its average. Because of universality it is possible to obtain exact results. In general, 
the average of an observable is not given by Random Matrix Theory. There are many 
examples of this type of universal fiuctuations ranging from atomic physics to quantum 
field theory (a recent comprehensive review was written by Guhr, Miiller-Groeling and 
Weidenmiiller^^). Most examples are related to fluctuations of eigenvalues. Typical 
examples arc nuclear spectra^^, acoustic spectra^^, resonances in resonance cavitics^'^, 
^-matrix fluctuations^^' and universal conductance fluctuations'^^. In these lectures 
we will discuss correlations in the bulk of Dirac spectra and the microscopic spectral 
density. The second type of application of Random Matrix Theory is as a schematic 
model of disorder. In this way one obtains qualitative results which may be helpful 
in understanding some physical phenomena. There are numerous examples in this 
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category. We only mention the Anderson model of localization^^, neural networks^^, 
the Gross- Witten model of QCD^"^ and quantum gravity^^ . In these lectures we will 
discuss chiral random matrix models at nonzero temperature and chemical potential. In 
particular, we will review recent work by Stephanov^^ on the quenched approximation 
at nonzero chemical potential. 

In the first lecture we will review some general properties of Dirac spectra including 
the Banks-Casher formula. From the zeros of the partition function we will show that 
there is an intimate relation between chiral symmetry breaking and correlations of 
Dirac eigenvalues. Starting from Leutwyler-Smilga sum-rules the microscopic spectral 
density will be introduced. We will discuss the statistical analysis of quantum spectra. 
It will be argued that spectral correlations of 'complex' systems are given by Random 
Matrix Theory. We will end the first lecture with the introduction of chiral Random 
Matrix Theory. 

In the second lecture we will compare the chiral random matrix model with QCD 
and discuss some of its properties. We will review recent results showing that the 
microscopic spectral density and eigenvalue correlations near zero virtuality are strongly 
universal. Lattice QCD results for the microscopic spectral density and correlations in 
the bulk of the spectrum will be discussed in detail. We will end the second lecture 
with a review of chiral Random Matrix Theory at nonzero chemical potential. Novel 
features of spectral universality in nonhermitean matrices will be discussed. 



THE DIRAC SPECTRUM 
Introduction 

The Euclidean QCD partition function is given by 

^QCD(m, 6) = J dAdet (^D + m)e-^^^/'^+^^^ (2) 

where ^yD — ^,j,{d^ + iA^) is the anti-Hermitean Dirac operator and Sym is the Yang- 
Mills action. The integral over field configurations includes a sum over all topological 
sectors with topological charge v. Each sector is weighted by exp(i^z/). Phenomenolog- 
ically the value of the vacuum 6'-angle is consistent with zero. We use the convention 
that the Euclidean gamma matrices are Hermitean with {'y^n^v) = 25^1,. The inte- 
gral is over all gauge field configurations, and for definiteness, we assume a lattice 
regularization of the partition function. 

Our main object of interest is the spectrum of the Dirac operator. The eigenvalues 
Afe are defined by 

-iD(t)k = i\k(t>k- (3) 

The spectral density is given by 

p(A) = $:5(A-A,). (4) 

k 

Correlations of the eigenvalues can be expressed in terms of the two-point correlation 
function 

P2(A,A') = (p(AKA')), (5) 
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where (■ • ■) denotes averaging with respect to the QCD partition function (^). The 
connected two-point correlation function is obtained by subtraction of the product of 
the average spectral densities 

p,(A,A')=P2(A,A')-(p(A))(p(A')). (6) 

Because of the Ua{^) symmetry 

{75,7^} = 0, (7) 

the eigenvalues occur in pairs ±A or are zero. The eigenfunctions are given by 0^ 
and 750fc, respectively. If •y^cpk = then necessarily A^ = 0. This happens for a 

solution of the Dirac operator in the field of an instanton. In a sector with topological 
charge u the Dirac operator has u exact zero modes with positive chirality. In order 
to represent the low energy sector of the Dirac operator for field configurations with 
topological charge u, it is natural to choose a chiral basis with n right-handed states 
and m = n + u left-handed states. Then the Dirac matrix has the block structure 

(lyt (8) 

where W is an n x m matrix. For m = 2 and n = 1, one can easily convince oneself 
that the Dirac matrix has exactly one zero eigenvalue. We leave it as an exercise to the 
reader to show that in general the Dirac matrix has \m — n\ zero eigenvalues. 

In terms of the eigenvalues of the Dirac operator, the QCD partition function can 
be rewritten as 

ZqcdK^) = 5:e-^n^/' / dAYlY[iXl + m})e-'--'/^ (9) 

u f f k 

where dA denotes the integral over field configurations with topological charge z/, 
and 11/ is the product over Nj flavors with mass irif. The partition function in the 
sector of topological charge v is obtained by Fourier inversion 

1 /■27r 

Zu{m) = — / dee-'^'Z^cT^im, 6). (10) 



27r 7o 

The fluctuations of the eigenvalues of the QCD Dirac operator are induced by the 
fluctuations of the gauge fields. Formally, one can think of integrating out all gauge 
field configurations for fixed values of the Dirac eigenvalues. The transformation of 
integration variables from the fields. A, to the eigenvalues, A^, leads to a nontrivial 
" Jacobian" . Universality in Dirac spectra has its origin in this " Jacobian" . 

The free Dirac spectrum can be obtained immediately from the square of the Dirac 
operator. For a box of volume Li x L2 x L3 x L4 one finds 

An^ = ( {--^Y + {--^Y + i^Y + (-^Y 1 ' , (11) 




where we have used periodic boundary conditions in the spatial directions and anti- 
periodic boundary conditions in the time direction. The spectral density is obtained 
by counting the total number of eigenvalues in a shell of radius AL/27r. The result is 

pfr-(A) ~ V\\ (12) 

For future reference, we note that in the generic case, when the sides of the hypercube 
are related by an irrational number, asymptotically, the eigenvalues are uncorrelated, 
i.e. 

P2(A,A') = (p(A))(p(A')). (13) 
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Figure 1. The free Dirac spectral density (dotted curve) and the spectral density of the 
Kogut-Susskind Dirac operator for a gauge field configuration generated with (i = 2.4 (histogram). 
Both spectral densities are on a 12"* lattice and are normalized to unit area. 



Two examples of Dirac spectra are shown in Fig. 1. The dotted curve represents 
the free Kogut-Susskind Dirac spectrum on a 12^ lattice with periodic boundary condi- 
tions in the spatial directions anti-periodic boundary conditions in the time direction. 
For an Ni x N2 x x A^4 lattice this spectrum is given by 

Here, Ui = 0, 1, ■ ■ ■ , [N,/2] [i = 1, 2, 3) and ^4 = 0, 1, ■ ■ ■ , [{N4 + l)/2]. The Kogut- 
Susskind Dirac spectrum for an SU{2) gauge field configuration with (3 = 2.4 on the 
same size lattice is shown by the histogram in the same figure (full curve). We clearly 
observe an accumulation of small eigenvalues. 



The Banks-Casher Relation 



The order parameter of the chiral phase transition, (ipip), is nonzero only below 
the critical temperature. As was shown by Banks and Casher^, (ipip) is directly related 
to the eigenvalue density of the QCD Dirac operator per unit four-volume 

E= 1(^^)1 = lim^^^. (15) 

It is elementary to derive this relation. The chiral condensate follows from the partition 
function (all quark masses are chosen equal), 

- 1 1 2m 

(iPij) = - lim yj;^dm log Z{m) = - lim - x'^+m^ ^' ^^^^ 

f k k 
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If we express the sum as an integral over the average spectral density, and take the 
thermodynamic limit before the chiral limit, so that many eigenvalues are less than 
m, we recover (|T3]). The order of the limits in (|T3]) is important. First we take the 
thermodynamic limit, next the chiral limit and, finally, the field theory limit. As can 
be observed from (|TBp the sign of {ipip) changes if m crosses the imaginary axis. 

An important consequence of the Bank-Casher formula ([T5|) is that the eigenvalues 
near zero virtuality are spaced as 



AA = l/p(0) = tt/SI/. (17) 

This should be contrasted with the eigenvalue spectrum of the non-interacting Dirac 
operator. Then eq. ([12|) results in an eigenvalue spacing equal to AA ~ l/V^^^. 
Clearly, the presence of gauge fields leads to a strong modification of the spectrum 
near zero virtuality. Strong interactions result in the coupling of many degrees of 
freedom leading to extended states and correlated eigenvalues. On the other hand, 
for uncorrelated eigenvalues, the eigenvalue distribution factorizes, and for A 7^ 0, we 
have p(A) ~ A^^-'"'"^ in the chiral limit, i.e. no breaking of chiral symmetry. One 
consequence of the interactions is level repulsion of neighboring eigenvalues. Therefore, 
the two smallest eigenvalues of the Dirac operator, ±Amin repel each other, and the 
Dirac spectrum will have a gap at A = with a width of the order of 1/Ty . 



Spectral Correlations and Zeros of the Partition Function 

The study of zeros of the partition function has been a fruitful tool in statistical 
mechanics'^''' In QCD, both zeros in the complex fugacity plane and the complex 
mass plane have been studied^^' Since the QCD partition function is a polynomial 
in m it can be factorized as (all quark masses are taken to be equal to m) 

Zq^CJi{m,d) = W{m-mk). (18) 

k 

Because configurations of opposite topological charge occur with the same probability, 
the coefficients of this polynomial are real, and the zeros occur in complex conjugate 
pairs. For an even number of flavors the zeros occur in pairs im^. In a sector with 
topological charge this is also the case for even Nj x v. The chiral condensate is 
given by 

1 11 
S(m) = - lim TT^dra log ^qcd("^, 0) = - lim • (^9) 

For an even number of flavors, S(m) is an odd function of m. In order to have a 
discontinuity at m = 0, the zeros in this region have to coalesce into a cut along the 
imaginary axis in the thermodynamic limit. 

In the hypothetical case that the eigenvalues of the Dirac operator do not fluctuate 
the zeros are located at = ±iXk- In the opposite case, of uncorrelated eigenvalues, 
the eigenvalue distribution factorizes and all zeros are located at iicx, where = (Af,). 
As a result, the chiral condensate does not show a discontinuity across the imaginary 
axis and is equal to zero. 

We hope to convince the reader that the presence of a discontinuity is intimately 
related to correlations of eigenvalues^^ . Let us study the effect of pair correlations for 
one flavor in the sector of zero topological charge. The fermion determinant can be 
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written as 



iUi^' + AD) = E ( ^ ) m''''-'H>^l • • • A^. (20) 



There are 



k 
21 



(2/-1)!! 



ways of selecting I pairs from • • • A^. The average of each pair of different eigenvalues 
is given by 

(A^,A^) = + C2, (21) 
where cr^ is the expectation value of A^ and C2 is the connected correlator 

^2 = (A^A^)-(A^)(A^), m^n. (22) 
This results in the partition function 



N 12 



*1 ' . , 

N \ k 



Z{m) = E E^'^"^-'^ \ k][2l\ - (23) 



fc=0 l=Q 



After interchanging the two sums, one can easily show that Z{m) can be expressed as 
a multiple of a Hermite polynomial 



Z{m) = (-f )^/^H^((a^ + m')/sJ-2C,). (24) 

In terms of the zeros of the Hermite polynomials, Zk-, the zeros of the partition function 
are located at 



yJ-2C2 - (25) 



rrik = Zk 

Asymptotically, the zeros of the Hermite polynomials are given by Zk ~ 7r/c/2\/ N (with 
integer ±A;). In order for the zeros to join into a cut in the thermodynamic limit, they 
have to be spaced as ~ This requires that 

C2 ~ (26) 

The density of zeros is then given by 

dk , , 

Nm. (27) 

dm 

We conclude that pair correlations are sufficient to generate a cut of Z{m) in the 
complex m-plane, but the chiral symmetry remains unbroken. Pair correlations alone 
cannot suppress the effect of the fermion determinant. 

Leutwyler-Smilga Sum Rules 

We have shown that pair-correlations are not sufficient to generate a discontinuity 
in the chiral condensate. In this subsection we start from the assumption that chiral 

symmetry is broken spontaneously, and look for consistency conditions that are imposed 
on the Dirac spectrum. As has been argued by Gasser and Leutwyler^° and Leutwyler 
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and Smilga^, in the mesoscopic range (0), the mass dependence of the QCD partition 
function is given by (for simphcity, all quark masses have been taken equal) 

Z^^{m,e)^l df/e"^^^^^Trf/e-/-/_ ^28) 

JU<^G/H 

The integral is over the Goldstone manifold associated with chiral symmetry break- 
ing from G to H. For three or more colors with fundamental fermions G/H = 
SU{Nf) X SU{Nf)/ SU{Nf). The finite volume partition function in the sector of 
topological charge u follows by Fourier inversion according to ([T0|). The partition func- 
tion for = is thus given by (128|) with the integration over SU{Nf) replaced by an 
integral over U{Nf). The case of Nf = 1 is particularly simple. Then only a U{1) 
integration remains, and the partition function is given by^ ^^£o("^) ~ loi^VTi). Its 
zeros are regularly spaced along the imaginary axis in the complex m-plane, and, in 
the thermodynamic limit, they coalesce into a cut. 

The Leutwyler-Smilga sum-rules are obtained by expanding the partition function 
Zy{m) in powers of m before and after averaging over the gauge field configurations and 
equating the coefficients. This corresponds to an expansion in powers of m of both the 
QCD partition function (0) and the finite volume partition function ( ]28|) in the sector 
of topological charge u. As an example, we consider the coefficients of in the sector 
with z/ = 0. This results in the sum-rule 

where the prime indicates that the sum is restricted to nonzero positive eigenvalues. 

The next order sum rules are obtained by equating the coefficients of order m^. 
They can be combined into 

We conclude that chiral symmetry breaking leads to correlations of the inverse eigen- 
values. However, if one performs an analysis similar to the one in previous section, it 
can be shown easily that pair correlations given by (|30[) do not result in a cut in the 
complex m-plane. Apparently, chiral symmetry breaking requires a subtle interplay of 
all types of correlations. 

For two colors with fundamental fermions or for adjoint fermions the pattern of 
chiral symmetry breaking is different. Sum rules for the inverse eigenvalues can be 
derived along the same lines. The general expression for the simplest sum-rule can be 
summarized as^^' 



(31) 



XI' Ai\u\ + id\miG/H) + l)/Nfy 

where dim{G/H) is the number of generators of the coset manifold in (p^). 

The Leutwyler-Smilga sum-rules can be expressed as an integral over the average 



spectral density and spectral correlation functions. For the sum rule (|29|) this results 



m 

1/2S2 J A2 4Nf' ^ ' 

If we introduce the microscopic variable 

u = A1/S, (33) 
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this integral can be rewritten as 



1 , , u ., du 1 , , 

The thermodynamic hmit of the combination that enters in the integrand, 

will be called the microscopic spectral densit'if. This limit exists if chiral symmetry 
is broken. Our conjecture is that ps{u) is a universal function that only depends on 
the global symmetries of the QCD partition function. Because of universality it can 
be derived from the simplest theory with the global symmetries of the QCD partition 
function. Such theory is a chiral Random Matrix Theory which will be introduced later 
in these lectures. 

We emphasize again that the ?7a(1) symmetry of the QCD Dirac spectrum leads 
to two different types of eigenvalue correlations: spectral correlations in the bulk of 
the spectrum and spectral correlations near zero virtuality. The simplest example of 
correlations of the latter type is the microscopic spectral density defined in (|35|). 

We close this subsection with two unrelated side remarks. First, the QCD Dirac 
operator is only determined up to a constant matrix. We can exploit this freedom 
to obtain a Dirac operator that is maximally symmetric. For example, the Wilson 
lattice QCD Dirac operator, , is neither Hermitean nor anti-Hermitean, but 
is Hermitean. 

Second, the QCD partition function can be expanded in powers of before or 
after averaging over the gauge field configurations. In the latter case one obtains sum 
rules for the inverse zeros of the partition function. As an example we quote. 



(36) 



u=0 

where we have averaged over field configurations with zero topological charge. 



SPECTRAL CORRELATIONS IN COMPLEX SYSTEMS 

Statistical Analysis of Spectra 

Spectra for a wide range of complex quantum systems have been studied both 
experimentally and numerically (a excellent recent review has been given by Guhr, 
Miiller-Groeling and Weidenmiiller^^). One basic observation is that the scale of varia- 
tions of the average spectral density and the scale of the spectral fluctuations separate. 
This allows us to unfold the spectrum, i.e. we rescale the spectrum in units of the local 
average level spacing. Specifically, the unfolded spectrum is given by 

Ar = (p(A'))rfA', (37) 

J —CO 

with unfolded spectral density 

Punf(A) = E'^(A-Ar)- (38) 

k 

The fluctuations of the unfolded spectrum can be measured by suitable statistics. 
We will consider the nearest neighbor spacing distribution, P{S), and moments of the 



9 



number of levels in an interval containing n levels on average. In particular, we will 
consider the number variance, S2(n), and the first two cumulants, 71 (n) and 72 (^^)- 
Another useful statistic is the A3(n)-statistic introduced by Dyson and Mehta"^^. It 
is related to T,2{n) via a smoothening kernel. The advantage of this statistic is that 
its fluctuations as a function of n are greatly reduced. Both E2(n) and A^ln) can be 
obtained from the pair correlation function defined as 

y2(A, A') = -(Punf(A)punf (A')) + (Punf(A))(punf (A')). (39) 

Analytical expressions for the above statistics can be obtained for the eigenvalues 
of the invariant random matrix ensembles. They are defined as ensembles of Hermitean 
matrices with Gaussian independently distributed matrix elements, i.e. with probability 
distribution given by 

P{H) ~ e-"^^^'^. (40) 

Depending on the anti-unitary symmetry, the matrix elements are real, complex or 
quaternion real. They are called the Gaussian Orthogonal Ensemble (GOE), the Gaus- 
sian Unitary Ensemble (GUE) and the Gaussian Symplectic Ensemble (GSE), respec- 
tively. Each ensemble is characterized by its Dyson index (3 which is defined as the 
number of independent variables per matrix element. For the GOE, GUE and the GSE 
we thus have P — 1, 2 and 4, respectively. 

Independent of the value of /3, the average spectral density is the semicircle, 

N 



(p(A)) = -V2^. (41) 

71 

Analytical results for all spectral correlation functions have been derived for each 
of the three ensembles^^ via the orthogonal polynomial method. We only quote the 
most important results. The nearest neighbor spacing distribution, which is known 
exactly in terms of a power series, is well approximated by 

PiS)^ S'cxpi-aaS^), (42) 

where a^g is a constant of order one. The asymptotic behaviour of the pair correlation 
function is given by^^ 

1^2(A,A0~^g^^ for p = 2, (44) 

cos27r(A-A0 1 + (7r/2) sin27r(A - AQ 
^^(-^'^^ 4(A-A0 + An^iX - \'Y ^ " ^^^^ 

The 1/(A — A')^ tail of the pair correlation function results in a logarithmic dependence 
of the asymptotic behavior of ^2(71) and A3(n), 

E2(n) ~ (2/77^/3) logn and A3(n) ~ /3E2(n)/2. (46) 

Characteristic features of random matrix correlations are level repulsion at short dis- 
tances and a strong suppression of fluctuations at large distances. 

For uncorrelated eigenvalues the level repulsion is absent and one finds 

P(5)=exp(-5), (47) 

and 

S2(n) = n and A3(n) = n/15. (48) 
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Spectral Universality 



The main conclusion of numerous studies of eigenvalue spectra of complex sys- 
tems is that spectral correlations of classically chaotic systems arc given by RMT^^. 
As illustration of this so called Bohigas conjecture, we mention three examples from 
completely different fields. The first example is the nuclear data ensemble in which 
the above statistics are evaluated by superimposing level spectra of many different 
nuclei^^. The second example concerns correlations of acoustic resonances in irregular 
quartz blocks^^. In both cases the statistics that were considered are, within experi- 
mental accuracy, in complete agreement with the GOE statistics. The third example 
pertains to the zeros of Riemann's zeta function. Extensive numerical calculations^^ 
have shown that asymptotically, for large imaginary part, the correlations between the 
zeros are given by the GUE. 

The Gaussian random matrix ensembles introduced above can be obtained^^ from 
two assumptions: i) The probability distribution is invariant under unitary transforma- 
tions; ii) The matrix elements are statistically independent. If the invariance assump- 
tion is dropped it can be shown with the theory of free random variables"'^ that the 
average spectral density is still given by a semicircle if the variance of the probability 
distribution is finite. For example, if the matrix elements are distributed according to a 
rectangular distribution, the average spectral density is a semicircle. On the other hand, 
if the independence assumption is released the average spectral density is typically not 
a semicircle. For example, this is the case if the quadratic potential in the probability 
distribution is replaced by a more complicated polynomial potential V{H). Using the 
supersymmctric method for Random Matrix Theory, it was shown by Hackenbroich 
and Weidenmiiller^^ that the same supersymmetric nonlinear cr-model is obtained for a 
wide range of potentials V(H). This implies that spectral correlations of the unfolded 
eigenvalues are independent of the potential. Remarkably, this result could be proved 
for all three values of the Dyson index. 

Several examples have been considered where both the invariance assumption and 
the independence assumption are relaxed. We mention H — > H + A, where A is an 
arbitrary fixed matrix, and the probability distribution of H is given by a polynomial 
V{H). It was shown by P. Zinn-Justin^^ that also in this case the spectral correla- 
tions are given by the invariant random matrix ensembles. For a Gaussian probability 
distribution the proof was given by Brezin and Hikami^'^. 

The domain of universality has been extended in the direction of real physical 
systems by means of the Gaussian embedded ensembles^^' The simplest example 
is the ensemble of matrix elements of n-particle Slater determinants of a two-body 
operator with random two-particle matrix elements. It can be shown analytically that 
the average spectral density is a Gaussian. However, according to substantial numerical 
evidence, the spectral correlations are in complete agreement with the invariant random 
matrix ensembles^^. 

A large number of examples have been found that fall into one of universality 
classes of the invariant random matrix ensembles. This calls out for a more general 
approach. Naturally, one thinks in terms of the renormalization group. This approach 
was pioneered by Brezin and Zinn-Justin^^. The idea is to integrate out rows and 
columns of a random matrix and to show that the Gaussian ensembles are a stable 
fixed point. This was made more explicit in a paper by Higuchi et al.^'^. However, 
much more work is required to arrive at a natural proof of spectral universality. 

Although the above mentioned universality studies provide support for the validity 
of the Bohigas conjecture, the ultimate goal is to derive it directly from the underlying 
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classical dynamics. An important first step in this direction was made by Berry^^. He 
showed that the asymptotics of the two-point correlation function is related to sum- 
rules for isolated classical trajectories. Another interesting approach was introduced by 
Andreev et al.^^ who were able to obtain a supersymmetric nonlinear sigma model from 
spectral averaging. In this context we also mention the work of Altland and Zirnbauer^^ 
who showed that the kicked rotor can be mapped onto a supersymmetric sigma model. 



CHIRAL RANDOM MATRIX THEORY 

Introduction of the Model 

In this section we will introduce an instanton liquid®^' inspired chiral RMT for 
the QCD partition function. In the spirit of the invariant random matrix ensembles 
we construct a model for the Dirac operator with the global symmetries of the QCD 
partition function as input, but otherwise Gaussian random matrix elements. The 
chRMT that obeys these conditions is defined by^' 

Nf 

Z^^,,(mi,---,m;v,) = / Wndetp + m;)e-^^'^^(^^^), (49) 



where 



I 



and is a 77. X m matrix with v = |77 — 777| and N = n + m. As is the case in QCD, we 
assume that u does not exceed \/N, so that, to a good approximation, n = N/2. The 
parameter is identified as the dimensionless volume of space time. The potential V 
is defined by 

V{<P) = Et^'- (51) 
k>i ^ 

The simplest case is the Gaussian case when V{(f)) = Below it will be shown 

that the microscopic spectral density is independent of the higher order terms in this 
potential provided that the average spectral density near zero remains nonzero. The 
matrix elements of W are either real {(3 = 1, chiral Gaussian Orthogonal Ensemble 
(chGOE)), complex {(3 = 2, chiral Gaussian Unitary Ensemble (chGUE)), or quaternion 
real {(3 = 4, chiral Gaussian Symplectic Ensemble (chGSE)). This partition function is 
invariant under 

W U^WV (52) 

where the n x n matrix U and the m x m matrix V are orthogonal matrices for (3 = 1, 
unitary matrices for (3 = 2, and symplectic matrices for (3 = 4. This invariance makes 
it possible to express the partition function in terms of eigenvalues of W defined by 

W = U^AV. (53) 

Here, A is a diagonal matrix with diagonal matrix elements > 0. In terms of the 
eigenvalues the partition function (|4^) is given by 

z^^,,(7n„ ■ ■ ■ , 777^,) = / dx\A{xi) f n xr^^'^^'-' n m>^i + 

k f f,k 

(54) 
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where the Vandermonde determinant is defined by 



A{Xl) = l[{Xl-X^). (55) 

k<l 

This model reproduces the following symmetries of the QCD partition function: 

• The Ua{'^) symmetry. All nonzero eigenvalues of the random matrix Dirac oper- 
ator occur in pairs ±A or are zero. 

• The topological structure of the QCD partition function. The Dirac matrix has 
exactly |z/| = |n — m| zero eigenvalues. This identifies u as the topological sector 
of the model. 

• The flavor symmetry is the same as in QCD. For /3 = 2 it is SU{Nf) x SU{Nf), 
for /3 = 1 it is SU{2Nf) and for /3 = 4 it is SU{Nf). 

• The chiral symmetry is broken spontaneously with chiral condensate given by 

E = lim 7rp(0)/iV. (56) 

AT— »oo 

[N is interpreted as the (dimensionless) volume of space time.) The symme- 
try breaking pattern is^^ SU{Nf) x SU{Nf)/SU{Nf), SU{2Nf)/Sp{Nf) and 
SU{Nf)/0{Nf) for /? = 2, 1 and 4, respectively, the same as in QCD^^ 

• The anti-unitary symmetries. For three or more colors with fundamental fermions 
the Dirac operator has no anti-unitary symmetries, and generically, the matrix 
elements of the Dirac operator are complex. The matrix elements W^i of the 
corresponding random matrix ensemble are chosen arbitrary complex as well (/3 = 
2). For Nc — 2, the Dirac operator in the fundamental representation satisfies 

[CT2K,t'yD] = 0, (57) 

where C is the charge conjugation matrix and K is the complex conjugation 

operator. Because, {CT2KY = 1, the matrix elements of the Dirac operator can 
always be chosen real, and the corresponding random matrix ensemble is defined 
with real matrix elements {13 = 1). For two or more colors with gauge fields in the 
adjoint representation the anti-unitary symmetry of the Dirac operator is given 

by 

[CK, i^D] = 0. (58) 

Because {CKy = —1, it is possible to rearrange the matrix elements of the Dirac 
operator into real quaternions. The matrix elements Wki of the corresponding 
random matrix ensemble are chosen quaternion real (/? = 4). 

Together with the invariant random matrix ensembles, the chiral ensembles are 
part of a larger classification scheme. Apart from the random matrix ensembles dis- 
cussed in this review, this classification also includes random matrix models for disor- 
dered super-conductors^^. As pointed out by Zirnbauer^^, all known universality classes 
of Hermitean random matrices are tangent to the large classes of symmetric spaces 
in the classification given by Cartan. There is a one-to-one correspondence between 
this classification and the classification of the large families of Riemannian symmetric 
superspaces^^. 
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Calculation of the Microscopic Spectral Density 

The joint eigenvalue distribution of the nonzero eigenvalues for zero masses follows 
immediately from the partition function For Nf flavors and topological charge v 
the result for arbitrary potential is given by^^ 

p,{K ■ ■ ■ , A.) = C,^n n I A^. - A,T n Xf^^'^^'-' exp f-^ V V{Xl)\ , (59) 



where Cp^n is a normalization constant. For (3 = 2 the average spectral density and the 
spectral correlation functions can be derived from (|59|) with the help of the orthogonal 
polynomial method^^. The orthogonal polynomials in the variable x = are deflned 
by 

/ da;x''e-"^(")p„"(x)P^(x) = KS^n- (60) 
Jo 

where 

a = Nf + \u\ for p = 2. (61) 

In the Gaussian case, ^(A) = S^A, the associated polynomials are the generalized 
Laguerre polynomials. That is why this ensemble is also known as the Laguerre 
ensemble^^' The Vandermonde determinant can be rewritten as A(xj) = detfci[a;jr']. 
By the addition of linear combinations of rows this determinant can be expressed in 
terms of the orthogonal polynomials (§0]), A{xi) ~ detki[Pi"Li{xk)]- By multiplying the 
two determinants one obtains 

P/3(xi, ■ ■ ■ , x„) = ^ det K'^ixk, xi), (62) 
n! 

where the kernel is deflned by 

71—1 1 

^"(^'1/) = E T^P."(x)P,"(y)(a;y)'^/2e-"(^(-)+^(^))/l (63) 
The average spectral density is obtained by integrating over X2, ■ ■ ■ , Xn, 

/n 
Y[dxkpix,X2,- ■ ■ ,Xn) = K"-{x,x). (64) 

k=2 

In terms of the original variables, A = y/x, the spectral density is given by 

p(A) =2A7r"(A2,A2). (65) 
In the Gaussian case the spectral density is thus given by 

p(A) = 2A2"+iE2»+2n"+ie-'^^'^' J2 ^Ll{\^T?n)Ll{\^Y?n), (66) 

k=o "fc 

where the are the generalized Laguerre polynomials. With the help of the Christoffel- 
Darboux formula the sum can be expressed into the rfi^ order Laguerre polynomial and 
its derivative (with L'^{x) = — L$^t^(a;)), 
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Here, we introduced the microscopic variable 

z = A^SA = 2nEA, (68) 

and used the exphcit expression for the normahzation of the Laguerre polynomials, 
hi = r(a + k + l)/k\. The microscopic limit is defined by 

and can be evaluated with the help of the asymptotic relation 

where Ja is the ordinary Bessel function of degree a. In order to take the this limit 
most conveniently, we substitute the recursion relations 

Ll{z)=Ll_,{z)+Ll-\z). (71) 
This results in the microscopic spectral density 

ps{z) = '-{Jl{z) - J,Uz)J,_,{z)). (72) 

From the asymptotic relation for the Bessel function 

/ 2 \ -^/^ vr TT 

Jv{z) ~ ( — ) cos(2; — —V — —) for z ^ oo (73) 

we find that 



TizJ ' 2 4- 



lim psix) = —. (74) 

X— >oo TT 

The result for the average spectral density follows from the asymptotic properties 
of the Laguerre polynomials. It given by the semicircular distribution 

p(A) = (nSV7r)y^4/S2-A2. (75) 

Notice that the microscopic limit of the average spectral density coincides with the 
asymptotic limit ( |7^ of the microscopic spectral density. 

The two-point correlation function is obtained by integrating the joint spectral 
density over all eigenvalues except two. The microscopic limit can again be expressed 
in terms of Bessel functions^^. 

The spectral density and the two-point correlation function were also derived 
within the framework of the supersymmetric method of Random Matrix Theory 

The calculation of the average spectral density and the spectral correlations func- 
tions for (3 = 1 and /3 = 4 is much more complicated. However, with the help of skew- 
orthogonal polynomials^^' exact analytical results for finite N can be obtained as 
well. 

The microscopic spectral density for SU{2) with fundamental fermions {(3 = 1) is 
given by^^ 

wJ2a{z)J2a~l{w) - Z J2a-l{z) J2a{w) 

^ {zw)^Hz^-w^) ' ^ ' 
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where a is the combination 

« = A^/ - ^ + Y- (77) 

The microscopic spectral density in the symmetry class with (3 = 4: was first calculated 
by Nagao and Forrester'''^. It is given by 

Psiz) = 2z'^ / duu^ / dv [J4a_i(2uf z) J4a(2n2;) — v J^a-iC^uz) J^aC^uv z)] (78) 
JO Jo 

with 4a = Nf + 2\u\ + 1. 

The spectral correlations in the bulk of the spectrum are given by the invariant 
random matrix ensemble with the same value of (3. For (3 = 2 this was already shown 
three decades ago by Fox and Kahn^^. For (3 = 1 and (3 = 4 this was only proved 
recently'^^. 

Duality between Flavor and Topology 

As one can observe from the joint eigenvalue distribution, for (3 = 2 the dependence 
on Nf and v enters only through the combination Nf+\v\. This allows for the possibility 
of trading topology for fiavors. In this section we will work out this duality for the 
finite volume partition function. This relation completes the proof of the conjectured 
expression^^ for the finite volume partition function for different quark masses and 
topological charge z/. 

For (3 = 2 the partition function ([5^ ) obeys the relation 

— — ~ ZNf+u,o{mi,---,mNf,0,---,0), (79) 

11/ 

where the argument of the last factor has i' zeros. 

As an example, the simplest nontrivial identity of this type is given by 

Zi,i(m) ~mZ2,o(m,0). (80) 

Let us prove this identity without relying on Random Matrix Theory. According to the 
definition ( pHj) we have 

Zi,i{m) ~ J c/^e^^e'"^^™^^ (81) 

and 

^2o(m,0)~/ dUe™''^'^^'''\ (82) 

Ju&(2) 

where M is a diagonal matrix with diagonal elements m and 0. The integral over U 
can be performed by diagonalizing U according to U = Uie^'^'^Ui^ , and choosing Ui 
and (pk as new integration variables. The Jacobian of this transformation is 

J^A^e'^"). (83) 

The integral over Ui can be performed using the Itzykson-Zuber formula. This results 
in 

-i<t>2\2 



^2,0 



e 



■m[cos(pi — cos</)2j 



gmyScos^i ^mV'£cos(p2-j (84) 
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Both terms in the last factor result in the same contribution to the integral. Let us 
consider only the first term ~ exp(ml^S cos0i). Then the integral over 02 has to be 
defined as a principal value integral. If we use the identity 

(cos 01 - COS 02) V tan((0i + 02)/2) 

the 02-integral of the term proportional to sin0i gives zero because of the principal 
value prescription. The term proportional to cos0i trivially results in Zi i. We leave 
it as an exercise to the reader to generalize this proof to arbitrary Nj and v. 

The group integrals in finite volume partition function ( p8|) were evaluated by 
Leutwyler and Smilga^ for equal quark masses. An expression for different quark masses 
was obtained by Jackson et al7^. The expression could only be proved for z/ = 0. The 
above duality can be used to relate a partition function at arbitrary u to a, partition 
function at u = 0. This completes the proof of the conjectured expression for arbitrary 
topological charge. 

UNIVERSALITY IN CHIRAL RANDOM MATRIX THEORY 

In the chiral ensembles, two types of universality studies can be performed. First, 
the universality of correlations in the bulk of the spectrum. As discussed above, they 
are given by the invariant random matrix ensembles. Second, the universality of the 
microscopic spectral density and the eigenvalue correlations near zero virtuality. The 
aim of such studies is to show that microscopic correlations are stable against deforma- 
tions of the chiral ensemble away from the Gaussian probability distribution. Recently, 
a number of universality studies on microscopic correlations have appeared. They will 
be reviewed in this section. 

We wish to emphasize that all universality studies for the chiral ensembles have 
been performed for (3 = 2. The reason is that (3 = 1 and (3 = 4 are mathematically 
much more involved. It certainly would be of interest to extend such studies to these 
cases as well. 

In addition to the analytical studies to be discussed below the universality of the 
microscopic spectral density also follows from numerical studies of models with the 
symmetries of the QCD partition function. In particular, we mention strong support 
in favor of universality from a different branch of physics, namely from the theory of 
universal conductance fluctuations. In that context, the microscopic spectral density of 
the eigenvalues of the transmission matrix was calculated for the Hofstadter^^ model, 
and, to a high degree of accuracy, it agrees with the random matrix prediction^^. Other 
studies deal with a class random matrix models with matrix elements with a diverging 
variance. Also in this case the microscopic spectral density is given by the universal 
expressions^^ (|T6D and (|T^). 

The conclusion that emerges from all numerical and analytical work on modified 
chiral random matrix models is that the microscopic spectral density and the correla- 
tions near zero virtuality exhibit a strong universality that is comparable to the stability 
of microscopic correlations in the bulk of the spectrum. 

Of course, QCD is much richer than chiral Random Matrix Theory. One question 
that should be asked is at what scale (in virtuality) QCD spectral correlations deviate 
from RMT. This question has been studied by means of instanton liquid simulations. 
Indeed, at macroscopic scales, it was found that the number variance shows a linear 
dependence instead of the logarithmic dependence observed at microscopic scales^^. 
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More work is needed to determine the point where the crossover between these two 
regimes takes place. 



Invariant Deformations of the Gaussian Random Matrix Ensembles 

In a first class of universality studies one considers probability distributions that 
maintain unitary invariance. In this case the joint probability distribution is given by 
(|59D. In this section we consider the simplest nontrivial case with only ai and a2 different 
from zero and present the proof of Akemann, Damgaard, Magnea and Nishigaki^^' for 
this case. We wish to point out that the general case only leads to minor complications. 

This case was first studied by Brezin, Hikami and Zee^^. They showed that the 
microscopic spectral density is independent of 02. A general proof valid for arbitrary 
potential was given by Akemann, Damgaard, Magnea and Nishigaki^^' The essence 
of the proof is a remarkable generalization of the identity for the Laguerre polynomials, 

lim L„(-) = Jo(2v^) , (86) 

to orthogonal polynomials determined by an arbitrary potential V . This relation was 
proved by deriving a differential equation from the continuum limit of the recursion 
relation for orthogonal polynomials. This proof has been extended to the microscopic 
correlation functions of all chiral ensembles in a recent work by Akemann, Damgaard, 
Magnea and Nishigaki^^. 

In the normalization Pfc(O) = 1, the orthogonal polynomials (0) satisfy the recur- 
sion relation 

xPu{x) = -n [Pk+i{x) - Pkix) - QkiPkix) - Pk^iix))] , (87) 

where 

Tk = — — and Qk = . (88) 

Pk+i rkhk-i 

Here, the coefficient of in Pk{x) is denoted by pk- The fractions and the normal- 
izations hk can be determined from the relations 



/ dxnV\x)e-''^^''^Pkix)Pk{x) = - dx— e-'^^^''^Pkix)Pk\x\ 
Jo Jo dx L 



1 



(89) 

PCO POO d V 

/ dxnV'{x)e-'^^^''^xPkix)Pk{x) - hk - 2khk = - / dx— e-"^(^)a;Pfc(x)Pfe(x 
^0 Jo dx L 



= 
(90) 



For the potential 



V{x) = aix + ^x^, (91) 



these relations reduce to 

naihk + na2rkhk{l + Qk) = 1, 

(92) 

naiVkhkil + Qk) + na2rl (hk+i + hk{l + Qkf + hk-iQl'j - hk - 2khk = 0. 

(93) 
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In order to proceed we take the continuum limit of these relations, i.e. k ^ oo and 
n ^ cxD at fixed k/n = t. With 

r{t) = r{k/n) + 0{l/n), 

h(t) = -hikin) + 0(l/n2), (94) 
n 

to leading order in 1/n, these recursion relations can be rewritten as 

ai + 2a2r(t) = j^, (95) 
air(t) + 3a2r^{t) = t. (96) 

The functions r{t) and h{t) are given by the solution of these equations. However, 
we do not need the explicit solution. A more useful property is that they satisfy a 
differential equation that does not depend on oi and 

r'{t)h{t)-2h'{t)r{t) = h\t). (97) 

Remarkably, as was shown by Akemann, Damgaard, Magnea and Nishigaki^^' this 
relation is valid for any polynomial potential. This relation is the essential ingredient 
in reducing the continuum limit of the recursion relation ( pTf ) to a Bessel equation. 

To leading order in 1/n the r.h.s. of the recursion relation (|87D is zero. We have 
to collect the terms of subleading order. In the continuum limit we may write 

rk-i = r{t) - -r'{t), 
n 

hk-i = -h{t) i^h'{t), 

Pk-i{x) = P{t, x) - lj^P{t, x), etc.. (98) 
This results in the differential equation 

„^.P(M) = -M<)|^|P(M). (99) 

We observe that the continuum limit of the recursion relation exists if we take at the 
same time the microscopic limit with fixed n^/x. If we introduce the new variables 

V = nX = n^fx, 
u{t) = = / ^= 100 

the differential equation reduces to the Bessel equation 

-^u^Piu, v"^) + v'^Piu, v"^) = 0. (101) 
u du du 

Here, and below we omit the argument t of u{t). The general solution of this Bessel 
equation is given by 

P(m, t;^) = AJq{uv) + BNq{uv) (102) 
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The constants are determined by the boundary condition that the orthogonal polyno- 
mials are normalized as Pn{0) = 1. This results in A = 1 and B = 0. 

The average spectral density follows from ( pB| ) and (p^). For a = (we omit the 
superscript a) we obtain 



n-l 



p(A) = 2Ae-'^^(^ ) -rPkiX')Pk{X'). (103) 



k=0 



hk 



With the help of the Christoffel-Darboux formula the sum over the polynomials can be 
expressed in P„_i(A^) and its derivative. This results in 

p(A) = 2Ae-"^(^^)^ {PLiiX')Pn{X') - P:(A2)Pn-i(A^)) , (104) 

where the prime denotes differentiation with respect to A^. The leading order terms 
contributing to the continuum limit of this equation cancel. The subleading terms 
follow from the third equation in (|98|) . In terms of the new variables we have, 

P„_i(A2) ~ P„(A2) - -^Piu, \') = P„(A^) - -^M^v), (105) 

nat nat 



PniX') = ^Pn{\') - n'^^Muv). (106) 



Taking both the continuum limit and the microscopic limit, eq. ( |104|) reduces to 
/AN 2^i^) ( It/ . d d ^ . . 1 d ^ . . d ^ , ,\ 

= (^-j^J',{uv)Muv) - J'^{uv)Uuv) + J',{uv)J',{uv)j . (107) 

All polynomials that enter in the precursor to this equation are of the n'th order. 
Therefore, we can put t = 1. If we also use the relations J'q{z) = —Ji{z) and J[ = 
— Ji{z)/z + Jo{z), and the definition of u (see eq. ( [1001) ), we obtain 

p(A) = ^Ar2V(l) iJo{uil)v)Jo{uil)v) + Ji(u(l)i;) Ji(n(l)t;)) . (108) 

For A — >■ cx) the microscopic spectral density should coincide with p(0). This allows us 
to identify m(1) as 

u{l) = (109) 

For the microscopic spectral density defined by (with z = NHX = 7rp(0)A, notice that 
N = 2n) 



we obtain 



P5(^) = |U'W + ^iW), (in: 



in agreement with the result (^) derived for the Gaussian chiral ensemble. This result 
is independent of the specific shape of the random matrix potential. This proves its 
universality. 
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Noninvariant Deformations of the Chiral Gaussian Random Matrix Ensem- 
bles 



In a second class of universality studies one considers deformations of the Gaus- 
sian random matrix ensemble that violate unitary invariance. In particular, one has 
considered the case where the matrix W in (^) is replaced by 

W^W + A, (112) 

whereas the probability distribution of W is Gaussian. Because of the unitary invari- 
ance, the matrix A can always be chosen diagonal. The simplest case with A = ttT 
times the identity was considered by Jackson et al}^ . This model provides a schematic 
model of the chiral phase transition. In this section we will give a detailed derivation of 
the average spectral density. The aim is to shown that the spectrum changes dramati- 
cally with variations of the temperature parameter. Nevertheless, it could be shown^^ 
that the microscopic spectral density is temperature independent in the broken phase. 
We will show that for large matrices, the average resolvent defined by 

G{z) = ^i:rl , ^ \ (113) 



2n \z + ie — T> 
obeys the cubic equation^^' (the parameter S = 1 in (^)) 

- 2zG'^ + G{z'^ - vr^r^ + i) _ ^ = o. (114) 

Here, V is the ensemble of random matrices 



with probability distribution given by 

P{W) = exp(-nTrWW"^). (116) 

The average over P(W) is denoted by the brackets (■■■). Because the operator (|115|) 
has only a finite support, it is possible to expand the resolvent in a geometric series in 
l/{z — K) for z sufficiently large. Here, K is the matrix 

One finds by inspection that G{z) satisfies 



where Q is the matrix 



z-H 

It should be clear that Q is block diagonal with the block structure 



S = {-^)- (119) 
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where In is the n x n identity matrix. Therefore, we find that G{z) = g. The average 
over W can be carried out immediately to give 

This yields the following matrix equation for g and h: 



z 

-nT 





2 



which leads to the two independent equations 

zg — TiTh = 1 + g 
zh - TrTg = gh . (123) 

Elimination of h yields the equation 



z-g 

Evidently, it can be rewritten as a cubic equation for g. By taking the trace of Q one 



zg-^^ = l + g'. (124) 

z- g 



obtains the announced cubic equation ( |1 14| ) for G{z). For completeness we mention 
that this derivation can be rewritten''^ in terms of the so called blues function^°. 
The average spectral density given by 

p(A) = --lmG{z = A) (125) 

TT 

is a semicircle at ttT = and splits into two arcs at ttT = 1. For the spectral density at 
zero one obtains p(0) = — vr^T^, and therefore chiral symmetry is broken for vrT < 1, 
and is restored above this temperature. In spite of this drastic change in average spectral 
density, it could be shown^^ with the help of a supersymmetric formulation of Random 
Matrix Theory that the microscopic spectral density does not depend on T. 

The super-symmetric method in the first paper by Jackson et al}^ is not easily 
generalizable to higher order correlation functions. A natural way to proceed is to 
employ the super-symmetric method introduced by Guhr^^. In the case of /3 = 2 this 
method results in an analytical expression for the kernel determining all correlation 
functions. This approach was followed in two papers, one by Guhr and Wettig^'' and one 
by Jackson et al}^. The latter authors studied microscopic correlation functions for A 
in (|112|) proportional to the identity, whereas Guhr and Wettig considered an arbitrary 
diagonal matrix A. It was shown that independent of the matrix A, the correlations 
are given by the Bessel kernel*^. Of course, a necessary condition on the matrix A is 
that chiral symmetry is broken. Guhr and Wettig also showed that correlations in the 
bulk of the spectrum are insensitive to A. 

The deformation W —>■ W + A with the probability distribution for W given by 
an arbitrary invariant potential has not yet been considered. We have no doubt that 
universality proofs along the lines of methods developed by P. Zinn-Justin^^ can be 
given. 



LATTICE QCD RESULTS 

Recently, the Dirac spectrum in lattice QCD received a good deal of attention. In 
particular, the connection between the topology of field configurations and the spectrum 
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of the Wilson Dirac operator has been studied in detail^''' Other studies are 

related to the connection between the Wilson Dirac spectrum and the localization 
properties of the eigenfunctions^^ . 

In this section we will focus ourselves on the spectral correlations of the lattice QCD 
Dirac operator. Both correlations in the bulk of the spectrum and the microscopic 
spectral density will be studied. Consistent with universality arguments presented 
above, we find that spectral correlations are in complete agreement with chiral Random 
Matrix Theory. 

Correlations in the Bulk of the Spectrum 

Recently, Kalkreuter^^ calculated all eigenvalues of the Nc — 2 lattice Dirac oper- 
ator both for Kogut-Susskind (KS) fermions and Wilson fermions for lattices as large 
as 12^. For the Kogut-Susskind Dirac operator, D^^ , we use the convention that it is 
anti-Hermitcan. Because of the Wilson-term, the Wilson Dirac operator, , is neither 
Hermitean nor anti-Hermitean. Its Hcrmiticity relation is given by = 75/^'^ 75. 

Therefore, the operator 75-D'^ is Hermitean. However, it does not anti-commute with 
75, and its eigenvalues do not occur in pairs iA^. 

In the case of SU{2), the anti-unitary symmetry of the Kogut-Susskind and the 
Wilson Dirac operator is given by^^' 

[L'^^,T2X] = 0, and [75L'^,75CXt2] =0. (126) 

Because 

(r2X)2 = -l, and (75CKT2)' = 1, (127) 

the matrix elements of the KS Dirac operator can be arranged into real quaternions, 
whereas the Wilson Dirac operator can be expressed into real matrix elements. There- 
fore, we expect that eigenvalue correlations in the bulk of the spectrum are described by 
the GSE and the GOE, respectively^^. The microscopic correlations for KS fermions are 
described by the chGSE. However, the microscopic correlations for Wilson fermions are 
not described by the chGOE but rather by the GOE. Because of the anti-unitary sym- 
metry, the eigenvalues of the KS Dirac operator are subject to the Kramers degeneracy, 
i.e. they arc double degenerate. 

In both cases, the Dirac matrix is tri-diagonalized by Cullum's and Willoughby's 
Lanczos procedure^'^ and diagonalizcd with a standard QL algorithm. This improved 
algorithm makes it possible to obtain all eigenvalues. This allows us to test the accuracy 
of the eigenvalues by means of sum-rules for the sum of the squares of the eigenvalues 
of the lattice Dirac operator. Typically, the numerical error in the sum rule is of order 
10-^ 

As an example, in Fig. 1 we show a histogram of the overall Dirac spectrum for 
KS fermions at /3 = 2.4. Results for the spectral correlations are shown in Figs. 2, 3 
and 4. The results for KS fermions are for 4 dynamical flavors with ma — 0.05 on a 
12^ lattice. The results for Wilson fermions were obtained for two dynamical flavors 
on a 8'^ X 12 lattice. For the values of f3 and k, we refer to the labels of the flgure. For 
P > 2.4, with our lattice parameters for KS fermions, the Dirac spectrum near zero 
virtuality develops a gap. Of course, this is an expected feature of the weak coupling 
domain. For small enough values of k, the Wilson Dirac spectrum shows a gap at A = 
as well. In the scaling domain the value of k is just above the critical value of k. A 
qualitative description of the Wilson Dirac spectrum can be obtained with a random 
matrix model with the structure of the Wilson Dirac operator^"^. 
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The eigenvalue spectrum is unfolded by fitting a second order polynomial to the 
integrated spectral density of a stretch of 500-1000 eigenvalues. The results for S2(n), 
A^^n) and P{S) in Fig. 2 show an impressive agreement with RMT predictions. The 
fluctuations in ^2(71) are as expected from RMT. The advantage of As-statistic is well 
illustrated by this figure. We also investigated^^ the n dependence of the first two 
cumulants of the number of levels in a stretch of length n. Results presented in Fig. 
4 show a perfect agreement with RMT. Spectra for different values of f3 have been 
analyzed as well. It is probably no surprise that random matrix correlations are found 
at stronger couplings. What is surprising, however, is that even in the weak-coupling 
domain (/3 = 2.8) the eigenvalue correlations are in complete agreement with Random 
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Figure 2. Spectral correlations of Dirac eigenvalues for Wilson fermions (left) and KS-fermions 
(right). Results are shown for the number variance, S2(n), the A3— statistic and the nearest 
neighbor spacing distribution, P{S). The solid, dashed and dotted curves represent the analytical 
result for the GOE, GUE and GSE, respectively. 



Matrix Theory. Finally, we have studied the stationarity of the ensemble by analyzing 
level sequences of about 200 eigenvalues (with relatively low statistics). No deviations 
from random matrix correlations were observed all over the spectrum, including the 
region near A = 0. This justifies the spectral averaging which results in the good 
statistics in Figs. 2 and 3. 
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In the case of three or more colors with fundamental fermions, both the Wilson and 
Kogut-Susskind Dirac operator do not possess any anti-unitary symmetries. Therefore, 
our conjecture is that in this case the spectral correlations in the bulk of the spectrum 
of both types of fermions can be described by the GUE. In the case of two fundamental 
colors the continuum theory and Wilson fermions are in the same universality class. 
It is an interesting question of how spectral correlations of KS fermions evolve in the 
approach to the continuum limit. Certainly, the Kramers degeneracy of the eigenvalues 
remains. However, since Kogut-Susskind fermions represent 4 degenerate flavors in the 
continuum limit, the Dirac eigenvalues should obtain an additional two-fold degeneracy. 
We are looking forward to more work in this direction. 




Figure 3. The number variance, T,2{n) and the first two cumulants, 71 (n) and 72 (rt) as a function of 
n for eigenvalues of the Wilson Dirac operator (left) and the Kogut-Susskind Dirac operator (right). 
The full and dotted curves represent the analytical result for the GOE and the GSE, respectively. 



The Microscopic Spectral Density 

The advantage of studying spectral correlations in the bulk of the spectrum is 
that one can perform spectral averages instead of ensemble averages requiring only a 
relatively small number of equilibrated configurations. This so called spectral ergodicity 
cannot be exploited in the study of the microscopic spectral density. In order to gather 
sufficient statistics for the microscopic spectral density of the lattice Dirac operator a 
large number of independent configurations is needed. One way to proceed is to generate 
instanton-liquid configurations which can be obtained much more cheaply than lattice 
QCD configurations. Results of such analysis^^ show that for Nc = 2 with fundamental 
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fermions the microscopic spectral density is given by the chGOE. For Nc = 3 it is given 
by the chGUE. One could argue that instanton-liquid configurations can be viewed 
as smoothened lattice QCD configurations. Roughening such configurations will only 
improve the agreement with Random Matrix Theory. 

Of course, the ultimate goal is to test the conjecture of microscopic universality 
for realistic lattice QCD configurations. In order to obtain a very large number of 
independent gauge field configurations one is necessarily restricted to relatively small 
lattices. The first study in this direction was reported recently^''' In this work, 
the quenched SU (2) Kogut-Susskind Dirac operator was diagonahzed for lattices with 
linear dimension of 4, 6, 8 and 10, and a total number of configurations of 9978, 9953, 
3896 and 1416, respectively. The results were compared with predictions from the 
chGSE. 
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Figure 4. The distribution of the smallest eigenvalue (left) and the microscopic spectral density 
(right) of the Kogut-Susskind Dirac operator for two colors and f3 = 2.0. Lattice results are 
represented by the histogram, and the analytical results for the chGSE are given by the dashed 



We only show results for the largest lattice. For more detailed results, including results 
for the two-point correlation function, we refer to the original work. In Fig. 4 we show 
the distribution of the smallest eigenvalue (left) and the microscopic spectral density 
(right). The lattice results are given by the full line. The dashed curve represents 
the random matrix results. The distribution of the smallest eigenvalue was derived by 
Forrester^'^ and is given by 

P(A„.in) = «y|(«A^i,)3/2/3/2(aA^in)e-^("'-'")', (128) 

where a = VE. The random matrix result for the microscopic spectral density is given 
in eq. (|78|). We emphasize that the theoretical curves have been obtained without any 
fitting of parameters. The input parameter, the chiral condensate, is derived from the 
same lattice calculations. The above simulations were performed at a relatively strong 
coupling of /? = 2. Recently, the same analysis^^ was performed for f3 = 2.2 and for 
(3 = 2.5 on a 16^ lattice. In both cases agreement with the random matrix predictions 
was found^*. 
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An alternative way to probe the Dirac spectrum is via the valence quark mass 
dependence of the chiral condensate^^ defined as 



S(m) 



1 

N 



dX{p{X)) 



2m 



\2 _|_ ^2 



(129) 



The average spectral density is obtained for a fixed sea quark mass. For masses well 
beyond the smallest eigenvalue, S(m) shows a plateau approaching the value of the 
chiral condensate S. In the mesoscopic range (|l]), we can introduce u = XmN and 
X = mNTj as new variables. Then the microscopic spectral density enters in S(m). For 
three fundamental colors the microscopic spectral density for (3 = 2 (eq. (ff^)) applies 
and the integral over A in ( |129| ) can be performed analytically. The result is given by^^, 



x{Iaix)Ka{x) + Ia+lix)Ka-l{x)), 



(130) 



where a = Nf + |z/| and la and Ka are modified Bessel functions. In Fig. 2 we plot 
this ratio as a function of x (the 'volume' V is equal to the total number of Dirac 
eigenvalues) for lattice data of two dynamical fiavors with mass ma = 0.01 and Nc = 3 
on a 16^ X 4 lattice. We observe that the lattice data for different values of (3 fall on a 
single curve. Moreover, in the mesoscopic range this curve coincides with the random 
matrix prediction for Nj = u = 0. 
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Figure 5. The valence quark mass dependence of the chiral condensate (rn) plotted as S^(to)/E 
versus mVYj. The dots and squares represent lattice results by the Columbia group -"^^ for values of /3 
as indicated in the label of the figure. 



Apparently, the zero modes are completely mixed with the much larger number of 
nonzero modes. For eigenvalues much smaller than the sea quark mass, one expects 
quenched (A^/ = 0) eigenvalue correlations. In the same figure the dashed curves 
represent results for the quark mass dependence of the chiral condensate (i.e. the mass 
dependence for equal valence and sea quark masses). In the sector of zero topological 
charge one finds^^' ^ 

Y.Hu) h(u) , , 
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and 



E»(«)__rfM_ _ 
— - u(iS(u) - iHu)) 

We observe that both expressions do not fit the data. Also notice that, according 
to Gockeler et a/.^°°, eq. ( |131[ ) describes the valence mass dependence of the chiral 
condensate for non-compact QED with quenched Kogut-Susskind fermions. However, 
we were not able to derive their result (no derivation is given in the paper). 



CHIRAL RANDOM MATRIX THEORY AT fi ^ 
Generalities 

At nonzero temperature T and chemical potential /i a schematic random matrix 
model of the QCD partition function is obtained by replacing the Dirac operator in 
(H) by^^' 101' 35 

iW + inT + fi\ 
^-[iW^ + tnT + fx )■ ^^^^^ 

Here, Qt = T C?>n (2n + l)7rl are the matrix elements of ijodo in a plane wave basis 
with anti-periodic boundary conditions in the time direction. Below, we will discuss 
a model with Qt absorbed in the random matrix and fi ^ 0. The aim of this model 
is to explore the effects of the non-Hermiticity of the Dirac operator. For example, 
the random matrix partition function (^) with the Dirac matrix ( |133| ) is well suited 
for the study of zeros of this partition function in the complex mass plane and in the 
complex chemical potential plane. Numerical results for the location of zeros could be 
explained analytically^"^. 

The term /i7o does not affect the anti-unitary symmetries of the Dirac operator. 
This is also the case in lattice QCD where the color matrices in the forward time 
direction are replaced by U —>■ e^U and in the backward time direction by —>■ e~^U'^ . 
For this reason the universality classes are the same as at zero chemical potential. 

The Dirac operator that will be discussed in this section is thus given by 



= ..„,,, (134) 




where the matrix elements of the n x n matrix W are either real (/5 = 1), complex 
(/3 = 2) or quaternion real (/3 = 4). For all three values of /3 the eigenvalues of T>{fi) 
are scattered in the complex plane. 

Since many standard random matrix methods rely on convergence properties based 
on the Hermiticity of the random matrix, direct application of most methods is not 
possible. The simplest way out is the Hermitizationi''^ of the problem, i.e we consider 
the Hermitean operator 



D''(z,z')- . (135) 




For example, the generating function in the supersymmetric method of Random Matrix 
Theoryio^' ig then given by^"''' 1°^' i°«' 

7(TT*\ / detjP'' {z + J, z* + .r) \ 
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The determinants can be rewritten as fermionic and bosonic integrals. Convergence is 
assured by the Hcrmiticity and by the infinitesimal increment k. The resolvent follows 
from the generating function by differentiation with respect to the source terms 



(137) 

,7=,7*=0 



Notice that, after averaging over the random matrix, the partition function depends in 
a non-trivial way on both z and z* . The spectral density is then given by 

Therefore, outside the domain of the eigenvalues G{z, z*) is a function of z only, whereas 
for z inside the domain of eigenvalues, G{z, z*) should depend on z* as well. 

Alternatively, one can use the replica trick^^*^' with generating function given by 

Z{z, z\ii) = {dei^WH^z, z*)). (139) 

The resolvent is then given by 

The idea is to perform the calculation for integer values of Nf and perform the limit 
Nf — * at the end of the calculation. Although the replica limit, A^^ ^ 0, fails in 
general^^^, it is expected to work for dei Dh{z,z*) because it is positive definite (or 
zero). Then the partition function is a smooth function of Nf. For other techniques 
addressing nonhermitean matrices we refer to the recent papers by Feinberg and Zee^*^^ 
and Nowak and co-workers^^^. One recent method that does not rely on the Hcrmitic- 
ity of the random matrices is the method of complex orthogonal polynomials'^^. This 
method was used by Fyodorov et al.^^ to calculate the number variance and the near- 
est neighbor spacing distribution in the regime of weakly nonhermitean matrices. As 
surprising new result, they found an 5"^/^ repulsion law. 

In the physically relevant case of QCD with three colors, the fermion determinant 
is complex for nonzero chemical potential. Its phase prevents the convergence of fully 
unquenched Monte-Carlo simulations (see Kogut et al}^^ for the latest progress in this 
direction). However, it is possible to perform quenched simulations. In such calculations 
it was found that the critical chemical potential //c ~ instead of a third of the 
nucleon mass^^^. This phenomenon was explained analytically by Stephanov^^ with the 
help of the above random matrix model. He could show that for small /i the eigenvalues 
are distributed along the imaginary axis in a band of width ~ ii^ leading to a critical 
chemical potential of Hc ~ V^- ^ detailed derivation of his solution will be given in 
the next section. As has been argued above, the quenched limit is necessarily obtained 
from a partition function in which the fermion determinant appears as 

lim \ <letV{^)ff, (141) 

instead of the same expression without the absolute value signs. The partition function 
with the absolute value of the determinant can be interpreted as a partition function of 
an equal number of fermions and conjugate fermions. The critical value of the chemical 
potential, equal to half the pion mass, is due to Goldstone bosons with a net baryon 
number consisting of a quark and conjugate anti-quark. The reason that the quenched 
limit does not correspond to the standard QCD partition function is closely related to 
the failure of the replica trick in the case of a determinant with a nontrivial phase. 
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The Stephanov Solution 



In this section we study the quenched limit of the partition function ( |139|) . We 
give a detailed derivation of the results originally obtained by Stephanov'^^. The deter- 
minants in (|139| ) can be written as Grassmann integrals. Since Grassmann integrals are 
always convergent the infinitesimal increment k can be put equal to zero. This results 
in the partition function 



Z{z,z*,ij,) = J VWVi)*Vil}V(t)*V(t)e-xY>[-nT?i:iWW^] 



X exp 



(142) 



where W is an arbitrary complex n x n matrix. The Gaussian integrals over W can be 
performed trivially, 



il^Ri'^^Ri^Lk - ^Lk'4^Ri(t>Ri(t>Lk ' A^Ri^Ri'^Lk + ^Lfe'/'R^Ri 



(143) 

The four-fermion terms can be written as the difference of two squares. Each square 
can be linearized by the Hubbard-Stratonovitch transformation according to 

exp(-AQ2) ^ / (i(jexp( -r - iQfj) ■ (144) 

J '-±/\ 

Using this, the fermionic integrals can be performed, and the partition function can be 
written as an integral over the complex 2Nf x 2Nf matrices, a and a, 

Z{z,z*,iJ,) = JvaVaexp[-nJ:'^Tr{a-C)k{a-C)k]det''(^aa-iJ,^). (145) 

In a Nf X Nf block structure notation, the matrices a, a, k and ( are given by 

_ ( a ih \ ^ _ ( ic^ \ 

where a, 6, c and d are arbitrary complex Nf x Nf matrices. The resolvent is obtained 
by differentiation with respect to z according to (|140|) with the averaged partition 
function given by ( |145| ). This results in 

G= -^Tr(2z + a + at). (147) 
INf 

In the thermodynamic limit the integrals in ( p.45|) can be evaluated by a saddle- 
point method. Notice that a variable and its complex conjugate have to be considered 
as independent integration variables. The saddle point equations are given by 

- ^ ^ d = k{d + C)k, (148) 



^ a = k{a + C)k. (149) 



ao — /i^ 
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In general the solution of the saddle point equations is not unique. However, a unique 
solution is obtained from the requirement that d*G{z) is positive definite. The saddle 
point equations have two obvious solutions 

a = a or a = kak, (150) 

for which the matrix equations (|148|) and ( |149|) coincide. 

It can be shown that the first possibility results in an unphysical solution. As 
usual in applications of the replica trick, we assume that the replica (or flavor) sym- 



metry remains unbroken. This implies that, at the saddle point, each block in ( |146| ) is 
diagonal. With (|150|) the solution of our saddle point equations is reduced to a 2 x 2 
matrix problem 

a = {akak - fi^){k(rk + (). (151) 
Consistent with the second solution in (|150|) we use the parametrization 

-=f (152) 



One solution of ( |151| ) can be written down immediately, namely b = b* = 0. Then 
the equations for a and d reduce to the same cubic equation. This solution results in 
a partition function which factorizes in a product of a z dependent part and in a z* 
dependent part. This resolvent is therefore an analytic function of z valid outside the 
domain of the eigenvalues. The cubic equation for the resolvent in this domain can be 
obtained from the cubic equation ( |1 14| ) by the replacement T —>■ in. 

Let us now focus on the solution inside the domain of eigenvalues with 6 7^ 0. 
From the off-diagonal elements of ( |151| ) we obtain the equations 



(a + d){d- z*) - (a^ - bb* - fx^) = 1, 

(a + rf)(a - z) - (rf2 - 66* -/i^) = 1. (153) 
The difference of these equation results in 

2{a + d){d-a) + {a + d){z- z*) = 0, (154) 

with solution given by 

a-d=^{z - z*). (155) 
From the sum of the two equations one obtains 

bb* = l + ha + d){z + z*)-ad-fi^. (156) 
The boundary of the domain of the eigenvalues is given by the set of points where this 



solution merges with the factorized solution, i.e. where 6 = 0. With (|155|) the condition 
6 = can be expressed as an equation for a, or by ( p.47| ) (G = a — z for diagonal a) 
an equation for the resolvent. On the boundary of the domain the equation for the 
resolvent is given by 

(2G + 2x + iy)x - (G + a; + i?/)(G + x) + 1 - /i^ = 0. (157) 
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For b = the Stephanov solution and the solution of the cubic equation merge into 
each other. We thus have a second order phase transition, and on the boundary, the 
second derivative of the free energy corresponding to ( |145D vanishes. 
From the diagonal elements of (|151D one obtains 

a = (a - z){a'^ - bb* - ji^) + bb*{a + rf), 

d={d- z*){d^ - bb* - I?) + bb*{a + d). (158) 



In the first equation we substitute the first equation of (|153|) in the first term of the 



r.h.s. and the expression ( |1 56| ) for bb* in the second term of the r.h.s.. This results in 
d-a + z + {a + d){^{z - z*){a - rf) - /i^ + zz*) = 0. (159) 
Together with ( |155| ) this results in an equation for a with solution 

+ ^Z/- (160) 



jj? — x"^ ' 2 

For the resolvent one obtains 



G = a-z = ^r^ — --x-\. (161) 



jj? — x"^ 



This results in the spectral density 




(162) 



inside the boundary given by 



x"^ 

2/2 = 4 - V + - — —. (163) 

yU^ — X'^ {jj,'^ — X'^Y 



The latter equation has obtained by the substitution of ( |161| ) in ( |157|) . Both ( |161 
and (|163|) were first derived by Stephanov'^^' The closed curves in Fig. 6 show 
the boundary of the domain of eigenvalues given by ( |163|) . Numerical results for the 
eigenvalues are represented by the dots in the same figure. 

For small /i, the width of the domain of eigenvalues is ~ /x^. This explains that 
the critical value of the chemical potential is given /ic ~ ^/m. This result explains the 
quenched results found in lattice QCD at nonzero chemical potential. 

Random Matrix Triality at Nonzero Chemical Potential 

In this section, we study the Dirac operator ( p.34| ) for all three values of (3. Both for 
(3 = 1 and /5 = 4 the fermion determinant, det(P(/i) + m), is real. This is obvious for 
(3 = 1. For (3 = 4: the reality follows from the identity q* = cr2Q'o"2 ^oi a quaternion real 
element g, and the invariance of a determinant under transposition. We thus conclude 
that quenching works for an even number of flavors. Consequently, chiral symmetry 
will be restored for arbitrarily small nonzero /i, whereas a condensate of a quark and 
a conjugate anti-quark develops. Indeed, this phenomenon has been observed in the 
strong coupling limit of lattice QCD with two colors^^^. 

In the quenched approximation, the spectral properties of the random matrix en- 
semble ( |134|) can be easily studied numerically by simply diagonalizing a set of matrices 
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Figure 6. Scatter plot of the real (x), and the imaginary parts (y) of the eigenvalues of the random 
matrix Dirac operator. The values of P and fj, are given in the labels of the figure. The full curve 
shows the analytical result for the boundary. 

with probability distribution (^). In Fig. 6 we show numerical results^^^ for the eigen- 
values of a few 100 x 100 matrices for fi = 0.15 and fi = 0.5. The dots represent 
the eigenvalues in the complex plane. The solid line is the analytical result^^ for the 
boundary of the eigenvalues which is given by the algebraic curve ( |163| ). This result 
was first derived by Stephanov for /? = 2 (see previous section). However, the method 
that was used can be extended^^^ to P = 1 and P = 4. Although the effective partition 
function is much more complicated, it can be shown without too much effort that the 
solutions of the saddle point equations are the same if the variance of the probability 
distribution is scaled as l/jS. In particular, the boundary of the domain of eigenvalues 
is the same in each of the three cases. However, as one observes from Fig. 6, for (3 = 1 
and (3 = 4 the spectral density deviates significantly from the saddle-point result. For 
(3 = 1 we find an accumulation of eigenvalues on the imaginary axis, whereas for (3 = 4 
we find a depletion of eigenvalues in this domain. This depletion can be understood as 
follows. For fi = all eigenvalues are doubly degenerate. This degeneracy is broken at 
fi ^ which produces the observed repulsion of the eigenvalues. 

The number of purely imaginary eigenvalues for /3 = 1 appears to scale as \/N. 
This explains that this effect is not visible in a leading order saddle point analysis. 
From a perturbative analysis of ( |139| ) one obtains a power series in 1/A^. Clearly, the 
\/N dependence requires a truly nonperturbative analysis of the partition function (|49|) 
with the Dirac operator ( p.34| ). Such a scaling behavior is typical for the regime of 
weak non-hermiticity first identified by Fyodorov et al}^^. Using the supersymmetric 
method for the generating function ( |136D the ^/N dependence was obtained analytically 
by Efetov^o^. 

A similar cut below a cloud of eigenvalues was found in instanton liquid simula- 
tions^^^ for A^c = 2 at /i 7^ and in a random matrix model of arbitrary real matrices^^"^. 
The depletion of the eigenvalues along the imaginary axis was observed earlier in lattice 
QCD simulations with staggered fermions^^°. Obviously, more work has to be done in 
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order to arrive at a complete characterization of universal features^^ in the spectrum 
of nonhermitean matrices. 



CONCLUSIONS 

We have argued that there is an intimate relation between correlations of Dirac 
eigenvalues and the breaking of chiral symmetry. In the chiral limit, the fermion deter- 
minant suppresses gauge field configurations with small Dirac eigenvalues. Correlations 
counteract this suppression, and are a necessary ingredient of chiral symmetry break- 
ing. From the study of eigenvalue correlations in strongly interacting systems, we have 
concluded that they are described naturally by a Random Matrix Theory with the 
global symmetries of the physical system. In QCD, this led to the introduction of chi- 
ral Random Matrix Theories. They provided us with an analytical understanding of 
the statistical properties of the eigenvalues on the scale of a typical level spacing. It 
could be shown analytically that the microscopic spectral density is strongly universal. 
These results constitute the foundation of the impressive agreement between lattice 
QCD and chiral Random Matrix Theory for the microscopic spectral density and for 
spectral correlations in the bulk of the spectrum. 

An extension of this model to nonzero chemical potential provided us with a com- 
plete analytical understanding of the failure of the quenched approximation observed 
in lattice QCD simulations at finite density. Some intriguing properties of previously 
obtained lattice QCD Dirac spectra and instanton liquid Dirac spectra at finite density 
could be explained as well. 
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